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PROCEDURES  FOR 

HEURISTIC  SaiEDlILING  UNDER  LIMITED  RESOURCES 
IN  ACTIVITY  NETWORKS 

Abstract 

This  is  the  documentation  of  a  program  that  permits  the  analyst 
to  choose  any  one  of  three  priority  rules  to  schedule  activities  under 
conditions  of  limited  resources.  The  program  is  written  in  FORTRAN. 

We  give  the  logic  of  the  programs,  a  flowchart  representation,  and  a 
printout  of  the  instructions  (upon  request).  Also,  a  solved  sample 
problem  is  included. 

This  program  is  available  at  a  nominal  fee  from  the  Graduate 
Program  in  Operations  Research,  N.  C.  State  University,  P.  0.  Box  5511, 
Raleigh,  NC  27650, 


PROCEDURES  FOR  HEURISTIC  SCHEDULING  UNDER 


LIMITED  RESOURCES  IN  ACTIVITY  NETWORKS 

S.  E.  Elmaghraby  and  Z,  M.  Naman 

This  package  consists  of  two  programs: 

1.  SETUP:  To  setup  the  variables  in  a  form  readable 

by  HSP  (heuristic  scheculing  procedure) 

2,  HSP:  To  schedule  the  activities  prepared  by  SETUP. 

I,  SETUP 

This  program  is  to  set  up  the  variables  needed  for  the  schedul¬ 
ing  program  HSP,  namely  the  network,  activities*  duration,  resource 
consumption  and  resources  availabilities.  The  following  are  its  options 
GN:  GN=1  The  network  is  randomly  generated 
GN=2  The  network  is  read  in 

GRAND:  GRAND=1  The  network  parameters  are  randomly  generated. 
GRAND=2  The  network  parameters  are  read  in 
CONV:  C0NV=1  Convert  the  network  to  A-on-N. 

C0NV=2  The  network  is  A-on-N- 
Any  combination  of  the  above  options  could  be  used.  For 
example. 


GN= 1  ,  GRAND=2  ,  CONV= 1 

means:  To  generate  the  network  randomly; 

The  network  parameters  are  read  in;  and 

The  network  is  A-on-A,  so  convert  it  to  A-on-N. 
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Input  Sequence 

1.  Options:  GN,  GRAND,  CONV 

2.  Number  of  nodes,  number  of  arcs 

3.  ^  the  network  is  read  THEN 

IF  it  is  A-on-A  THEN  give  a  list  of  pairs  (Node  i. 
Node  j)  to  represent  arcs  in  the  network 
ELSE  [The  network  is  A-on-N]  give  the  upper  triangular 
matrix  without  diagonal  (the  adjacency  matrix) 
which  is  a  (0,1)  matrix  given  row  by  row. 

ELSE  Generate  the  network  randomly; 

4.  Read  the  number  of  resources 

5.  the  network  parameters  are  randomly  generated  THEN  give 
the  following  information: 


1. 

L3: 

Maximum 

desired 

activity  duration 

2. 

N3: 

Maximum 

desired 

resource  requirement 

3. 

F  >  1: 

Resource  allocation  factor 

[To  get  the  resource  allocation  we  find  the  average 
resource  used,  then  we  multiply  it  by  F. 

Average  resource  =  Total  resource/critical  path.] 
ELSE  [The  parameters  are  read  in]  give  the  following 
information: 

1.  Duration  of  each  activity 

2.  Resource  requirements  to  be  given  as  FOR  each 
activity 


FOR  each  resource 


Rl;  Resource  required 

Dl:  Duration  of  this  request 

1:ND 

Resource  availability  for  each  resource 
END 

Example;  1,  1,  1  [option  card]  [generate  the 

network,  generate  the  parame¬ 
ters,  convert  to  A-on-N] 

20,  60  [network  specification  card] 

r20-nodes  and  60-arcs] 

3  [number  of  resources  card] 

[3-resources] 

15,  10,  1,5  [parameters  generated  informa¬ 
tion]  [L3=15,  N3=10,  F=l,5] 

Intermediate  Output : 

1.  Provided  or  generated  information 

2.  Network  as  A-on-N  (arc  list) 

3.  Number  of  activities 

4.  Upper  bound  on  project  duration 

5.  Lower  bound  on  project  duration 

6.  Activity  data 

7.  Resource  availability  levels. 

Program  Limitations: 

1.  Maximum  number  of  resources  =  4 

2.  Maximum  number  of  activities  =  99 

3.  Maximum  sum  of  durations  =  350 

To  increase  the  above  variables,  the  program  DIMENSION 
must  be  changed. 
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4,  To  add  the  network  generator,  do  the  following: 

a.  Replace  the  comment  card  which  indicates  the 
position  of  the  network  generator  by  calling 
card. 

b.  Take  out  all  the  **c"  at  Column  1  in  the  next 
Do- loop . 

c.  Add  the  arrays  NE  and  NS  in  the  DIMENSION 
statement . 

5.  For  the  random  number  generator,  please  look  at  the 
note  given  in  "NETWORK  GENERATION"  program  descrip¬ 
tion,  OR  Report  No.  179,  June  1981,  p.  4. 

II.  HSP 

The  project  scheduling  problem  is  to  assign  start  times  to  all 
the  activities  in  A,  Each  of  these  start  times  must  respect  the  pre¬ 
cedence  and  resource  constraints  described  above.  Another  assumption 
made  in  this  study  is  that  of  non -preemptive  scheduling:  Once  started, 
an  activity  must  continue  without  interruption  for  its  duration.  As 
well  as  respecting  all  constraints,  schedulers  attempt  to  assign  start 
times  so  that  some  objective  is  achieved. 

The  objective  of  the  scheduling  algorithm  presented  below  is  to 
minimize  the  project  duration.  Define  the  set,  A(i),  to  be  the  set  of 
activities  immediately  succeeding  i.  By  L  we  denote  the  set  of  activi¬ 
ties  having  no  successors.  Activity  i  is  in  L  if  A(i)  is  the  empty 

set.  Project  duration  is  then  obtained  as  the  max  (t(i)  +  d(i)  -  Ij  . 

ieL 

Necessarily,  the  maximization  of  total  resource  utilization  is  accom¬ 
plished  by  minimizing  project  duration.  Moreover,  minimizing  project 


duration  also  assures  the  minimization  of  delay  penalties  incurred 
per  period  beyond  the  (resource  free)  critical  path  length.  Finally, 
the  chosen  objective  is  applicable  to  a  variety  of  situations.  For 
example,  a  capital  intensive  project,  such  as  a  large  construction 
project,  might  require  substantial  financing  to  secure  the  necessary 
resources.  Interest  payments  on  the  borrowed  capital  cut  into  the 
constructor's  profit  until  the  project  is  completed  and  the  resources 
are  no  longer  needed.  Minimizing  project  duration  thus  minimizes 
such  interest  costs. 

Heuristics  Employed 

In  the  heuristic  scheduling  procedure,  an  adapted  binpacking 
rule  is  used  to  assign  the  activities  to  the  resource  units.  The 
first  fit  decreasing  (FFD)  binpacking  rule  is  one  of  the  most  popular 
scheduling  rules.  It  has  also  been  applied  with  some  success  to 
the  problem  of  memory  storage  allocation  by  some  computer  scientists. 
The  FFD  rule  specifies  the  packing  of  different  sized  objects  into 
containers  or  bins.  The  objects  are  first  ordered  by  decreasing  size. 
The  largest  object  is  then  assigned  to  the  first  container  in  which 
it  will  fit.  The  second  largest  object  is  then  assigned  to  the  first 
container  in  which  it  will  fit,  etc.  In  the  project  setting,  the 
objects  are  the  jobs  and  the  containers  are  the  resources.  There  are 
as  many  "containers”  as  there  are  different  resources,  with  each 
job  "occupying"  all  containers  simultaneously.  This  is  a 
multi -dimensional  binpacking  problem  which  we  wish  to  avoid  by 


jj 


adopting  a  simple  rule  for  "sizing"  the  jobs.  We  address  this 
"sizing"  problem  next. 

For  project  scheduling,  it  is  clear  that  the  limited  avail¬ 
ability  of  the  various  resources  may  cause  an  activity  to  be  delayed. 
Such  delays  can  result  in  project  completion  delay.  A  measure  of 
the  delay  in  project  completion  caused  by  delaying  an  activity 
should  serve  well  as  an  activity's  "size," 

The  original  idea  for  measuring  the  project  delay  caused  by 
delaying  activity  i  was  to  sum  the  durations  of  all  activities 
succeeding  activity  i  in  the  precedence  ordering  plus  the  duration 
of  activity  i  itself.  The  larger  the  priority  number  ("size")  of 
activity  i,  p(i),  the  longer  the  expected  project  delay  should  be, 
if  activity  i  is  delayed.  To  reduce  the  computational  burden,  we 
adopted  a  slightly  different  priority  number: 


d(i) 

if  A(i)  =  (j). 

d(i)  + 

S  PO). 

if  A(i)  ^  (j). 

jeA(i) 


These  numbers  can  be  assigned  in  one  pass  through  the  precedence 
ordering.  However,  the  recursion  allows  the  counting  of  some  dura¬ 
tions  more  than  once  in  determining  the  priority  number  for  certain 


activities . 


7 


The  two  other  priority  rankings  available  in  this  package 
are:  the  longest  remaining  path  (LRP)  heuristic  and  the  shortest 
imminent  activity  (SIA)  heuristic.  LRP  gives  a  higher  priority  to 
an  activity  with  a  longer  critical  subpath  (to  project  completion) 
than  a  competing  activity.  The  SIA  heuristic  mimics  the  shortest 
processing  time  rule  from  job  shop  scheduling.  Characteristic  of 
both  schemes,  and  a  number  of  other  heuristics,  is  the  neglect  of 
the  total  impact  on  project  completion  that  delaying  an  activity  can 
produce.  The  impact  of  delaying  an  activity  can  affect  all  its 
successors.  The  proposed  priority  ranking  is  an  attempt  to  reflect 
the  full  impact  of  delaying  an  activity  on  the  project  duration  by 
adopting  a  simple  surrogate  measure. 

To  sum  up,  the  three  heuristic  "sizing^^  priority  rules  are: 

CRl :  The  sum  of  durations  of  all  succeeding 
activities 

CR2:  The  longest  path  to  project  termination  (LRP) 

CR3:  The  shortest  imminent  (available)  activity  (SIA) 

The  scheduling  process  begins  by  forming  the  set,  CS,  of  all 
activities  where  predecessors  have  been  scheduled  prior  to  the  current 
scheduling  period.  On  the  first  iteration,  CS  will  consist  of  all 
initial  activities  (i  such  that  B(i)  =  $) .  In  general  CS  will  contain 
activities  that  are  unordered  with  respect  to  precedence.  The  adopted 
FFD  rule  assigns  activities  of  CS  to  available  resource  units.  The 
activities  are  the  items,  ranked  by  decreasing  ^*size”  (priority  number). 
The  ’’size’*  of  each  bin  in  a  particular  scheduling  period  will  be  the 
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number  of  available  units  of  that  resource  type.  The  activities  are 
checked  for  fit  starting  with  that  activity  whose  priority  number  is 
the  largest  in  CS.  If  activity  i  *'fits",  then  activity  i  is  scheduled 
to  begin  on  period  Tl.  After  updating  resource  availabilities,  the 
next  activity  is  checked  for  fit,  etc.,  until  all  activities  in  CS 
have  been  considered  for  scheduling  at  Tl,  where  Tl  is  the  current 
scheduling  period. 

The  algorithm  is  iterative  in  nature.  CS  and  Tl  are  updated 
after  each  iteration.  The  process  repeats  until  all  activities  are 
assigned.  The  updated  scheduling  period  is  determined  by  the  con^^letion 
times  of  the  activities  scheduled  during  the  previous  iteration.  It  is 
the  earliest  completion  time  of  any  of  the  activities  scheduled  in  the 
previous  iteration. 

Tl  =  Tl  ,  .  +  min  (M(i,  Tl  .  '  d(i)), 

new  old  ^  ^  ’  old'^  j  ^  y 


where  d(i)  is  the  duration  of  activity  i,  and 


1  if  activity  i  started  on 


“  if  i  not  started  on  Tl  ,  ,, 

old 


The  intention  in  skipping  over  the  intervening  periods  between 
and  is  to  avoid  wasting  time  attempting  to  schedule  activities 

when  few  resources  are  available  in  the  scheduling  period.  Upon 
completion  of  an  activity,  resource  units  are  "released**  and  other 
activities  can  perhaps  be  undertaken.  A  flowchart  of  HSP  logic  is 
shown  in  Fig,  1. 


An  Example 


/T2.3) 


(6,3) 


(1,3) 


(3,3) 


(d(i),  c^^^) 


(7,3) 


(4,3) 


(5,3) 


The  figure  above  is  a  graphical  representation  of  the  precedence 
relationship  among  the  seven  activities  of  an  example  project.  An  arc, 
(i,j)  implies  that  activity  j  cannot  be  started  until  activity  i  has 
been  completed.  For  example,  activities  two,  three,  and  four  cannot 
be  started  until  activity  one  has  been  completed.  Assume  that  activity 
durations,  d(i)  =  i,  for  i  =  1,2,  ...,  7.  Also,  one  resource  type  is 
required  in  quantities  of  three  units  per  day  by  each  activity 

=  3,  for  j  =  1,2,  d(i)).  Six  units  of  the  resource  are 

available  for  each  of  the  scheduling  periods  (r(t)  =  6,  for 
t  =  1,2,  ...,  T  where  T  is  an  upper  bound  on  project  duration) . 

First,  the  priority  numbers  are  assigned  in  one  backward  pass 
through  the  precedence  network.  Suppose  we  adopt  priority  rule  CRl , 
Activities  six  and  seven  are  the  terminal  activities  of  this  project 
(i.e.,  L  =  {6,7}).  Thus  their  priority  numbers  are  just  their  durations 
(p(6)  =  dC6)  =  6,  pC7)  =  d(7)  =  7) .  Priorities  for  the  immediate  pre¬ 
decessors  of  six  and  seven,  activities  two,  three,  and  five,  can  now 
be  computed.  Their  priorities  are  the  sum  of  activity  duration  and 
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the  priority  numbers  of  the  successors:  p(2)  =  dC2)  +  p(6)  =  2  +  6  =  8, 
P(3)  =  d(3)  +  p(6)  =3  +  6  =  9,  p(5)  =  d(5)  +  p(7)  =  5  +  7  =  12. 
Activity  four  has  priority  p(4),  given  by  the  sum  of  its  duration  and 
the  priority  number  of  its  successor,  activity  five.  Thus  p(4)  = 
d(4)  +  p(5)  =  4  +  12  =  16,  Finally,  activity  one  has  as  successors 
activities  two,  three,  and  four.  Thus  p(l)  =  d(l)  +  p(2)  ♦  p(3)  + 
p(4)  =1+8+9+  16  =34.  Now  that  activity  priorities  have  been 
computed,  the  scheduling  process  can  begin. 

On  the  first  scheduling  iteration,  CS  will  consist  of  activity 
one  since  it  is  the  lone  initial  activity  in  the  project.  Three  units 
of  the  resource  are  required  by  activity  one,  and  six  units  are  avail¬ 
able  on  day  one.  Thus  activity  one  is  scheduled  for  day  one. 

Day  two  is  the  updated  scheduling  period  on  the  second  iteration 
since  activity  one  is  completed  on  day  one.  CS  will  consist  of  activi¬ 
ties  two,  three,  and  four,  the  immediate  successors  of  activity  one. 

The  scheduling  priority  is  four  (p(4)  =  16),  three  (p(3)  =  9),  and 
two  (p(2)  =  8).  Activity  four  is  first  checked  for  fit.  Three  units 
of  the  resource  are  required  on  days  two  through  five,  and  six  units 
are  available.  Since  activity  one  is  completed  by  day  two,  activity 
four  is  scheduled  to  start  on  day  two.  Now  only  three  resource  units 
remain  available  on  days  two  through  five.  Next,  activity  three  is 
checked  for  fit.  Three  resource  units  are  required  by  activity  three 
on  days  two  through  four,  and  three  units  are  available.  Activity  one 
is  completed  by  day  two,  so  activity  three  is  scheduled  to  start  on 
day  two.  Thus  no  units  of  the  resource  are  available  on  days  two, 
three,  and  four.  Activity  two  is  checked  for  fit  next.  But,  with  no 


resources  available  and  three  units  required,  activity  two  cannot  be 
scheduled  on  day  two.  So  the  set  CS  and  scheduling  period  are  updated 
and  we  proceed  to  the  third  iteration. 

The  new  set  of  activities  to  be  scanned  will  include  activity 
two  that  was  not  scheduled  on  the  previous  iteration  and  activity  five 
since  its  predecessor  has  now  been  scheduled.  The  new  scheduling  period 
is  day  five  since,  of  the  two  activities  scheduled  on  the  previous 
iteration,  activity  three  finishes  sooner  (on  day  four)  releasing  some 
resources.  Activity  five  has  a  larger  priority  number  than  activity 
two.  Hence  activity  five  is  considered  for  scheduling  before  activity 
two.  However,  activity  four  is  still  in  progress  on  day  five  so  activ¬ 
ity  five  is  ineligible  for  scheduling  on  day  five.  Next,  activity  two 
is  considered  for  scheduling  on  day  five.  Three  units  of  the  resource 
are  required  by  activity  two  on  days  five  and  six.  Three  units  are 
available  on  day  five,  and  six  units  are  available  on  day  six.  No 
precedence  constraint  is  violated,  so  activity  two  is  scheduled  to 
begin  on  day  five.  With  that  assignment,  no  resource  units  remain 
available  on  day  five  and  three  units  remain  available  on  day  six. 

Both  activities  in  CS,  at  this  iteration,  have  been  checked.  We  pro¬ 
ceed  to  update  CS  and  the  scheduling  period  for  the  next  iteration. 

The  updated  CS  consists  of  activities  five  and  six  since  the  scheduling 
of  activity  two  releases  activity  six  for  consideration.  The  new 
scheduling  period  will  be  day  six  since  activity  four  is  completed  on 
day  five. 

On  the  fourth  iteration,  activity  five  has  priority  and  is 
scheduled  on  day  six.  Available  resource  units  are  updated  and 


activity  six  is  considered.  Sufficient  resources  are  available,  and 
no  precedence  constraint  is  violated  in  starting  activity  six  on  day 
seven.  Hence,  that  assignment  is  made. 

Finally,  activity  seven  is  considered  on  the  fifth  iteration. 
The  updated  scheduling  period  is  day  eleven  (given  by  the  completion 
of  activity  five) ,  No  resource  or  precedence  constraints  are  violated 
by  starting  activity  seven  on  day  eleven  so  that  assignment  is  made. 
The  resource  "skyline"  or  profile  chart  below  illustrates  the  example 
project's  schedule.  The  project  is  completed  on  day  seventeen.  This 
duration  coincides  favorably  with  the  project's  critical  path  length 
in  the  absence  of  resource  constraints  (seventeen  days) ,  and  is  there¬ 
fore  optimal . 


The  project  scheduling  algorithm  described  above  is  an  analog 


of  the  heuristic  algorithm  for  the  ALB  problem  called  the  "Ranked 
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Positional  Weight”  heuristic,  which  was  proposed  by  Hegelson  and  Birnie 
in  1961,  An  extension  of  this  heuristic  was  suggested  by  Mansoor  in 
1964. 


Advantages  of  Heuristic  Scheduling 

A  significant  contribution  of  heuristic  scheduling  might  be  that 
of  increasing  the  computational  efficiency  of  optimal  project  schedul¬ 
ing  procedures.  One  attempt  might  be  to  link  the  proposed  scheduling 
algorithm  (or  any  other  good  heuristic  algorithm)  with  a  branch  and 
bound  algorithm.  For  such  a  tandem  algorithm,  the  size  of  the  space 
to  be  searched  will  be  smaller  than  the  search  space  of  a  straight¬ 
forward  branch  and  bound  procedure. 

We  have  noted  that  a  lower  bound  on  project  duration  is  the 
critical  path  length.  Similarly,  we  obtain  a  lower  bound  on  the  start 
time  of  an  activity  by  computing  the  earliest  start  schedule  (ESS) . 

The  ESS  neglects  any  resource  constraints  and  starts  an  activity  as 
soon  as  all  precursory  activities  are  completed: 


ESS(i)- 


r 


< 


max 

jeB(i) 


(ESS(j)  +  d(j)). 


if  B(i)  4  ^  . 


1  if  B(i)  =  (|l. 


If  a  tight  upper  bound  on  optimal  start  time  is  available  for  all 
activities,  then  the  chance  of  finding  an  optimal  schedule  by  a 
branch  and  bound  scheme  is  greatly  improved.  A  good  heuristic 
algorithm  will  provide  such  a  bound  at  a  small  price. 


The  heuristic's  schedule  yields  a  feasible  project  duration,  T. 


The  late  start  schedule,  LSS^,  analogous  to  the  early  start  schedule 
described  above,  is  computed: 


LSS^(i) 


min  (LSS  (j)  -  d(i)  +  1),  if  ACi)  4  i 
jeA(i) 


T  -  dCi)  +  1,  if  A(i)  =  (|). 


The  optimal  start  timeCs)  of  activity  i  is  (are)  now  bounded: 

ESS(i)  LSS^(i) .  Activity  i  cannot  be  started  before  its 

predecessors  are  completed  (ESSCi))  and  an  optimal  schedule  cannot 

start  activity  i  after  its  late  start  time  with  project  duration  T 

since  T  <  T. 
opt  — 

Thus  the  search  tree  of  the  branch  and  bound  procedure  is  pruned. 
The  tightness  of  these  bounds  and  hence  the  reduction  of  the  search 
space  is  increased  as  the  heuristic  schedule  is  improved.  Therefore, 
an  efficient,  dependable  heuristic  scheduling  algorithm 
might  play  a  crucial  role  in  increasing  the  computational  reliability 
of  an  analytical  procedure  for  the  project  scheduling  problem. 

The  following  information  is  needed  to  be  passed  from 

SETUP: 

1 ,  Niimber  of  resources 

2,  Number  of  activities 

3,  Upper  bound  of  planning  horizon 

4,  Network  in  A-on-A 

5,  Availability  of  each  resource 

6,  Requirement  of  each  resource/activity/time  period 


7.  Duration  of  each  activity 

8.  Index  vector  K 


9.  Early  start 

schedule 

Options 

P7; 

P7=l 

Do  not  print  iteration  information. 

P7«2 

Print  iteration  information. 

SCH: 

SCH=1 

Scheduling  using  CRl 

SCH=2 

Scheduling  using  CRl  and  CR2 

SCH=3 

Scheduling  using  CRl  and  CR3 

SCH=4 

Scheduling  using  CRl,  CR2,  and  CR3 

For  example,  P7=2 
SCH=3 


imply:  Print  iteration  information. 

Use  CRl  and  CR3  as  priority  criteria. 

Input  Sequence 

SCH  [on  the  same  card] 

Output 

1.  For  each  scheduling  iteration  (optional) 

a.  Current  scanned  subset  of  activities  and  their 
priority  numbers 

b.  Activities  delayed  because  of  resource  limited 
availability 

c.  Activities  delayed  because  of  precedence 

d .  Activities  scheduled 

2.  Project  table 

For  each  activity  Print: 

Duration,  Priority  number,  HSP-start,  Early-start 

J 
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3 .  Resource  utilization 
For  each  resource  Print: 

a.  Time,  availability,  utilization,  leftover 

b.  Draw  the  skyline  for  the  resource  utilization 
vs.  time. 

4.  Adjusted  percent  resource  utilization,  for  all  resources 
5-  Lower  bound  on  duration 

6.  Completion  time  from  HSP 

7.  Upper  bound  of  duration 
Limitation 

Please  see  the  limitation  given  in  SETUP. 


(OPTION)  To  read  the 
options  for  this  RUN, 

Read  all  the  necessary 
information  depending  on 
the  option  that  is  given, 

P  =  Number  of  nodes  if 
A-on-N 

=  Number  of  arcs  if 
A-on-A 

K  =  An  array  to  tell  where 
each  resource  require¬ 
ment  is  started  for 
each  activity - 


OPTION _ 

READ,  GN,  GRAND,  CONV _ 

Print  the  option  used  in  this  run 

READ,  NARC,  NODE _ 

Sv  GN  = 


NETGEN 


K  =  1 


C0NV=1 


NUP  =  2*NARC  NUP  =  NODE* 


FOR  1=  1  to  NARC  READ,  V(J)  (N0DE-l)/2 


V(K)  =  NE(I) 
V(K+1)  =  NS(I) 
K  =  K  +  2 


J  =  1,  NUP  J  =  1,  NUP 


READ,  RI 
GRAND  =  2 
READ,  L3,  N3,  F 
^  CONVERT 


READ,  D(I),  I 

=  1.  P 

I  =  1  to  P 

FOR  RO  =  1,  RI 

READ,  RI,  D1 

FOR  DO  =  1,  D 

1 

R(DO+K(I),RO 

)  =  RI _ 

_ READ,  A(1,R0 

FOR  I  =  1,  P 
B2  =  B2  -»•  D(I) 

FOR  RO  =  1,  RI _ 

FOR  I  =  2,  B2 
A(I,RO)  *  A(1,R0) 
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2 .  (CONVERT)  To  convert 
A-on-A  to  A-on-N  and 
create  an  array  X  which  is 
the  one  dimensional  repre¬ 
sentation  of  the  upper 
triangle  w/o  diagonal 
matrix. 

X  =  1  if  activity 
I  precedes  J 
=  0  otherwise 

3.  (GEN)  To  generate  randomly 
the  network  parameters. 

D(I)  -  duration  of  each 

activity 

R  -  resource  requirement 
for  each  activity  at 
each  time  period  for 
each  resource. 

4.  (CPL)  To  calculate  the 
critical  path  of  the  network. 
The  early  start  schedule 
could  be  produced  by 

E(J)  =  V(J)  -  D(J)  +  1; 

1  J  P  where  V  is  the 
vector  generated  by  CPL. 


CONVERT 


CONV  =  2 


FOR  I  =  1,  P  -  1 


IXl  =  (I-l)*(P-I/2)  -  I 


FOR  J  =  I  +  1,  P 


IX  =  1X1  +  J 
X(IX)  =  0 


V(2*I)  4  V(2*J-1) 


X(IX)  =  1 


GEN 


GRAND=2 


FOR  I  =  1,  P 


RANDU  (RAND) 


D(I)  =  [L3*RANDfl] 


FOR  I  =  1  to  P 


FOR  RO  =  1  to  RI 


FOR  DO  =  1  to  D(I) 


RANDU  (RAND) 


R(DO-t-K(I),RO)  =  fN3*RND-«-l] 


CPL 


FOR 

13  =  1  to  P 

V(I3)  =  D(I3)  j 

FOR 

I  =  1  to  P  -  1 

1X1  =  CI-l)*P-I/2)  -  I 

FOR  J  =  I  +  1  to  P 

IX  =  1X1  +  J 

K  X(IX)  =  0 

V(J)  =  MAXtV(J),  V(I) 

KPATH  =  V(P) 
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1. 


2. 


3. 


(INI)  Read  option 
variables  P7,SCH.  Set  the 
matrix  L  =  matrix  A,  Set 
vector  X  to  zero. 

(PRNO)  Compute  the  priority 
number  by  CRl .  Store  it  in 
vector  E. 

(COMPET)  Get  a  set  of 
simultaneous  activities. 

Save  the  candidates  in 
vector  S. 

X  =  0  if  the  activity  has 
not  been  scheduled 
X  =  T1  the  time  the  activity 
is  scheduled 
LL  =  0  No  event  between 

activities  PO  and  J 
LL  =  1  Otherwise 


M  =  M  1 
S(M)  =»  PO 
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4.  (ORDER)  Sort  the  set  S 
according  to  their  priority 
number . 

5.  (SCH)  The  beginning  of 
scheduling  the  set  S.W-  will 
keep  track  of  those  activities 
that  are  scheduled  or  cannot 
be  scheduled  because  of 
precedence  constraints . 

T1  is  the  time  counter. 

6.  (PRECON)  To  check  if  there 

is  any  precedence  constraints . 
IF  there  is  none  THEN  check 
for  resource  constraints, 

ELSE  go  to  CHEK. 

7.  (ADQ)  To  check  if  the 
resources  requested  by 
activity  S(W)  is  available. 

IF  yes  THEN  schedule  that 
activity,  and  update  the 
resource  available ,  ELSE 
go  to  CHEK. 


ORDER _ 

FOR  I  =  1  to  M  -  1 _ 

FOR  J  =  1  to  M  -  I 
\  E(SCI))  >  E(S(I-hJ)) 
N  =  S(I) 

S(I)  =  S(I+J) 


- - - — - - 1 

SCH 

w  = 

1 

y\^  T1  = 

1 

ADQ 

PRECON 

BL  =  0  ,  B2  = 

0  ,  PO  =  scw) 

FOR  J  =  1,  PO  -  1 _ 

PREC  (LL) 

Nv _  LL  =  0 _ 

_ B1  =  B1  1 

Vy  X(J)  =  0 

X(J)  D(J)  >  T1 

_ I  B2  =  B2  1 

_ B2  =  B1 

-*■  CHEK 


FOR  DO  =  1  to  D(S(W)) _ 

FOR  RO  =  1  to  RI _ 

R(DO-t-K(SW),RO)  :<.L(Tl-*-DO-l,RO) 

_ I  -*•  CHEK _ 

X(SW)  =  T1  ,  V(J2)  =  S(W)  ,  J2  =  J2  +  1 

FOR  DO  =  1,  D(S(W)) _ 

FOR  RO  =  1,  RI _ 

L(T1+D0-1,R0)  =  L(T1+D0-1,R0)- 
_  _  R(IXH-K,RO) 
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8.  (CHEK)  To  check  if  all 
the  candidate  activities 
are  checked.  ^  they  are, 
THEN  check  if  any  activity 
has  not  been  scheduled. 

IF  there  are  any,  THEN 
order  the  already  scheduled 
activities  and  prepare  to  go 
to  GOBACK.  ELSE  go  to 
OUTPUT.  ELSE  go  to  PRECON. 

9.  (GOBACK)  To  prepare  the 
variables  for  going  back 
either  to  schedule  whatever 
activities  are  left  in  S  or 
to  find  a  new  competing  set 
of  activities.  Update  the 
time  clock. 


CHEK 


W  =  M 

/n 

w  =  W  +  1 

^  PRECON 

FOR  PO  =  2,  P 

X 

o 

II 

o 

_ A 

FLAG  =  1 

FLAG  =  1 

_ A 

-►  OUTPUT 

K 

J2  =  0 

/n 

FOR  I  =  1  to  J2  -  1 

FOR  J  =  I  +  1  to  J2  -  1 

v(n  >  vcj) 

_ A 

N  =  V(I) 

V(I)  =  V(J) 

V(J)  ^  N 

GOBACK 


CO  =  1 

FOR  I  =  1  to  M 

YN.  X(S(I))>  0 

Cl  =  Cl  +  1 

\  Cl  4  M 

T1  =  T1  +  1 

SCH 

K3  =  0 

FOR  I  =  1  to  M 

X(S(I))  =0  A 

K3  =  K3  +  1 

K3  4  1  A 

T1  =  MAX(T1,X(S(I))>T1  -  T1  + 

+  D(S(I))  DCSCD) 

COMPET 

(OUTPUT)  To  print  the 
calculated  schedule. 

PREC  -  A  block  used  to 
calculate  the  vector  image 
of  the  adjacency  matrix. 
Then  pass  LL. 

MAIN  -  The  main  program 
that  calls  the  subroutines 
SETUP  and  HSP. 


OUTPUT 


IXl  =  (J-l)*(P-J/2)  -  J  *  PO 

_ LL  =  XB(IXI) _ 

RETURN 


SETUP 
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